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Abstract 


This  semi-annual  report  is  divided  into  three  separate 
sections  the  first  summarizes  the  conclusions  reached  in 
the  experimental  flow  discharge  facility  concerning  the 
deposition  of  energy  into  the  internal  modes  of  product 
molecules  in  simple  combustion  reactions,  the  second  summarizes 
the  progress  in  the  theoretical  close  coupling  vibrational 
relaxation  cross  sections  appropriate  to  rarified  flow 
regimes  and  the  last  section  reviews  the  work  to  date  in 
the  molecular  beam  experiments  aimed  at  determining  the 
interaction  potential  for  plume  species. 
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Part  I:  INFRARED  CHEMILUMINESCENCE  FROM  OXYGEN  ATOM 

ATTACK  ON  ETHYLENE 


R.H.  Kummler,  Mazen  Malki,  and  E.R.  Fisher 


Part  I 


Abstract 

The  2.7u  emission  from  oxygen  atom  attack  on  ethylene 
has  been  studied  in  a  discharge- flow  tube  system  under  condi¬ 
tions  of  low  ethylene  concentrations  and  excess  oxygen  atoms 
which  favor  the  production  of  vibrationally  excited  OH^  at 
total  pressures  of  about  1  Torr.  The  reactions  responsible  for 
this  emission  have  been  found  to  be 

CHO  ♦  O  -*>  CO  ♦  OH*  (A) 

followed  by 

OH^  -*■  OH  ♦  hv  (7000A  to  2.7y) 


Another  looser  but  potentially  important  icacuuu  lung  uuuldcl 
times  which  may  contribute  to  the  OH^  is 


CH20  ♦  0  •*  CHO  ♦  OH 


♦ 


(B) 


14 


An  upper  limit  for  kg  has  been  determined  to  be  3  x  10 

The  major  loss  mechanism  for  CHO  in  this  system  must  be: 


CHO  ♦  02  -*•  (products)  (C) 

and  the  ratio  of  the  rate  constants  for  reactions  (A)  and  (C)  was 
found  to  be  k^/k^,  ■  1.7  +  0.2. 
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Introduction 

The  energy  released  in  an  exothermic  chemical  reaction 
can  take  several  forms  depending  on  the  reaction  path.  The 

chemical  energy  of  the  reactants  can  be  converted  into  trans¬ 
lational  energy  of  the  products  or  it  can  be  converted  into 

internal  energy  of  the  products  or  some  combination  of  both. 

The  internal  energy  can  be  distributed  among  the  electronic, 
vibrational,  or  rotational  modes  of  the  product  molecule  or 
molecules.  If  the  product  is  elevated  to  a  quantum  state  which 
has  an'  optically  allowed  transition  to  a  lower  state,  then  the 
resulting  light  emitted  permits  the  study  of  the  reaction  path 
leading  to  the  internal  energy.  The  observation  of  radiation 
associated  with  this  non-equilibrium  process  thus  provides  in¬ 
formation  regarding  the  fundamental  chemical  kinetics  of  the 
reaction  (e.g.,  Charters  and  Polanyi  (1960);  Polanyi,  1956, 

Polanyi,  1971;  Hushafar,  et  al,  1971).  It  also  provides  the 
opportunity  to  assess  the  presence  of  one  or  both  of  the  reactant 
partners  in  an  unknown  mixture  of  reactant  gases.  This  capabil¬ 
ity  is  of  considerable  interest  in  the  air  pollution  monitoring 
field  and  is  applicable  to  both  local  and  remote  monitoring  of 
exhaust  gases  (of  rockets  as  well  as  automobiles)  in  the  presence 
of  an  oxygen  atom  environment.  In  the  spectral  region  of  wave¬ 
lengths  from  ly  to  roughly  15  0r  20U»  the  emitted  light  generally 
characterizes  vibrational  transitions  which  can  result  when  the 
exothermic  reaction  produces  a  new  bond  in  a  heteronuclear  molecule. 

One  common  product  of  the  combustion  of  hydrocarbons  is 
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ethylene;  oxygen  atom  attack  on  ethylene  will  produce  a  sequence 
of  reactions  between  secondary  species  eventually  producing 
several  heteronuclear  molecules  with  new  bonds,  especially  OH 
and  H2O,  which  are  strong  radiators  in  the  near  infrared.  Thus, 
we  might  expect  to  observe  light  emitted  at  the  fundamental 

vibrational  frequency  of  OH  as  well  as  the  mode  of  H^O,  both 
near  2.7p.  One  of  the  early  results  of  this  work  was  the  ob¬ 
servation  of  emission  at  that  wavelength.  In  order  to  apply 
the  results  of  such  observations  to  the  sensing  of  the  gaseous 
environment  of  interest,  it  is  necessary  to  establish  the 
molecular  mechanism  leading  to  the  radiation  because  both  the 
absolute  concentrations  as  well  as  the  ratio  of  reacting 
molecules  in  the  field  experiment  may  be  considerably  different 
than  in  the  laboratory.  With  this  in  mind,  the  major  task  of 
this  work  was  to  establish  the  radiation  contributing  reactions 
and  the  major  loss  mechanisms  for  the  observed  infrared  radiation 
In  order  to  avoid  competing  reactions,  it  was  decided  that  low 
concentrations  of  ethylene  would  be  employed  thereby  avoiding 
reactions  dependent  upon  more  than  one  of  the  hydrocarbon  frag¬ 
ments  and  in  order  to  restrict  the  system  to  one  radiator,  vib- 
rationally  excited  OH  (denoted  hereafter  by  011^  without  regard 
to  the  details  of  the  excitation).  This  process  is  necessary 
initially  to  simplify  the  reaction  mechanism  to  a  tractable 
point.  The  second  major  objective  then  becomes  the  quantifi¬ 
cation  of  the  reaction  system  in  terms  of  rate  constants  or 
quantum  yields  (ratio  of  photons  emitted  per  molecule  reacted) 
for  the  system. 
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The  ethylene -oxygen  atom  system  was  chosen  for  this 
study  because  the  ground  state  chemistry  of  the  system  is 
well  known  from  the  work  of  Niki,  et  al  (1969),  Stuhl  and 
Niki  (1971),  Kanofsky,  et  al  (1972),  Atkinson  and  Cvetanovic 
(1972)  and  the  intensive  effort  by  the  National  Bureau  of 
Standards  culminating  in  the  very  accurate  establishment  of 
the  rate  constant  for  the  primary  step 

0  +  C2H4  -►  CHO  +  CH3 

by  Davis,  et  al  (1972).  Under  these  favorable  conditions,  the 
concentrations  of  the  major  species  can  be  accurately  modeled 
in  the  experimental  system,  and  the  rate  constants  for  the 
radiating  step  can  be  obtained  by  essentially  a  perturbation 
analysis.  This  system  was  also  chosen  because  of  previous 
work  by  Krieger,  Malki,  and  Kummler  (1972)  which  identified 
OH^  as  a  major  radiator  in  the  very  near  infrared  (7000-9000A) 
under  identical  experimental  conditions. 
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EXPERIMENTAL  PROCEDURE 


The  oxygen  atom  attack  on  ethylene  under  1  torr  of 
total  pressure  was  studied  using  a  discharge  flow  system 
shown  schematically  in  Figure  (1). 

The  flow  tube  reactor  is  a  4  ft.  long  stainless -steel 
tube  with  a  1  inch  I.D.  teflon  tube  liner  which  aids  in  pre¬ 
venting  wall  recombination  of  some  atomic  species.  It  is 
provided  with  six  ports  for  raii-  j.  viewing  through  KC1  windows 
(Harshaw  Chemical  Company) .  Two  additional  ports  were  used 
for  measuring  the  absolute  pressure  and  the  drop  in  pressure 
across  the  reactor.  Linear  average  velocities  from  25  -  35 
m/sec  (with  33  m/sec  being  a  typical  velocity)  were  achieved  when 
thp  pressure  in  the  reactor  was  about  1  torr  using  a  Heraeus 
Englehard  (E-225,  147  cfm)  mechanical  pump.  A  Baratron 
(MKS  type  77)  capacitance  manometer  was  employed  to  monitor 
the  reactor  pressure  to  better  than  ly  accuracy. 

The  detectors  used  for  this  work  were  the  PbS  detector 
(Santa  Barbara  Research  ATO  type  8471-1)  in  conjunction  with 
a  phase  sensitive  lock-in  amplifier  (PAR  HR8)  for  IR 
detection,  and  a  photomultiplier  detector  (RCA  C31025C) 
in  conjunction  with  an  electrometer  (Keithley  602)  for 
observation  of  visible  emission. 


MKS  Baratrcn 
pressure  meter 
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In  order  to  get  absolute  emission  values,  the  system 


was  calibrated  using  the  reaction  (Fontijn  et  al ,  1964) 

NO  +  0  ■*>  N02  +  hv  (2 .  7y) 

for  converting  the  arbitrary  lock-in  amplifier  readings  from 
y  volts  to  absolute  photons/cc  sec.  The  actinometric  standard  of 
Fontijn,  et  al  was  extended  into  the  2.7y  region  using  the  data 

of  Stair  and  Kennealy  (1967). 

The  oxygen  atoms  were  generated  by  passing  oxygen 

(Cryogenic  sales  99.6%)  pure  or  diluted  in  argon  (Cryogenic  sales 

991)  through  an  Evanscn  type  cavity  discharge  powered  by  a  2450 

MH7  Ravtherm  PCM  TOY?  diathermy  unit. 

* 

The  oxygen  atoms  were  passed  into  the  flow  tube  reactor 
through  a  closed  end  teflon  tube  with  a  radially  multi -perfo 
rated  exit  stream.  This  permitted  thorough  mixing  and  a  plug 
flow  pattern  at  the  point  of  mixing  due  to  the  turbulence  caused 
by  introducing  the  gas  through  a  "showerhead"  gas  inlet.  The 
plug  on  the  end  of  the  teflon  piece  served  also  to  prevent  the 
discharge  region  radiation  from  reaching  the  observation  ports 
through  scattering.  The  oxygen  atom  concentrations  were 
measured  by  using  the  N02  titration  technique  (Kaufman  1961)). 

It  was  found  by  using  very  small  contrations  of  ethylene 

(on  the  order  of  <  ly)  complex  hydrocarbon-hydrocarbon  interactions 
could  be  eliminated  significantly  simplifying  the  kinetics  and 
making  the  analysis  of  this  system  possible.  A  calibrated 
capillary  provided  the  needed  restriction  on  the  C2H4  flow  into 
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the  reactor. 


A  steady  state  flow  of  the  gas  appropriate  to  the  experi¬ 
ment  C 0 2  or  O2  +  Ar)  was  obtained.  Argon  served  to  keep  the 
pressure  in  the  reactor  constant  allowing  a  reasonably  constant 

linear  velocity  upon  adding  the  ground  state  gas. 

Thus,  Argon,  from  a  side  gas  inlet,  contributes  most  of  the 
pressure  in  the  reactor  (°1  terr)  while  th 


rt  rrr  At  m  A  r  f  a  +■  a  o  f  k»a  a 

v  L  a  w  mum  bu  ^  bi  vum 


(NO,  NO^ ,  C2H4)  and  the  discharged  stream  (80y  0^)  are  relatively 
small  fractions  of  the  total  pressure. 

The  phase  angle  on  the  lock-in  amplifier  changed  as  one 
moved  the  PbS  detector  downstream  and  as  the  pressure  changed 
in  the  flow  tube.  This  had  to  be  checked  for  all  experimental 
data  points.  The  reference  to  the  PAR  was  taken  from  the  pulse 
generator  of  the  PGM  10X2.  The  pulse  rate  of  the  discharge  was 
adjusted  so  that  a  single  pulse  of  atoms  traveled  down  the  flow 
tube.  This  typically  required  repetition  rates  of  SO  Hz  to  100 
Hz.  The  duty  cycle  for  most  measurements  was  501. 

A  run  that  involved  the  reaction  NO  +  0  -*•  N02  +  hv  accom¬ 
panied  every  set  of  data  to  put  the  lock-in  amplifier  readings 
on  an  absolute  basis.  A  titration  for  oxygen  atoms  was  also  done 
as  a  part  of  the  system  calibration.  The  PbS  detector  with  the 


PAR  IIR8  was  used  to 


monitor  th  £  intensity  peak  of  the  leactiuii 


N02  +  0  NO  +  02  as  a  function  of  N02  concentration,  a  sample 


titraticr  is  given  in  Figure  2. 

All  observations  of  the  0  +  C2H4  emission  were  made  using 
the  2.7y  filter  with  transmission  characteristics  shown  in 
Figure  3. 
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Partial  Pressure  of  N0o,  p 

%  L. 

Figure  2.  Titration  of  0  using  NC^  ,  observing  the  NO  +  0 
glow  with  the  PbS  detector  and  the  2  .  7p  filter  removed. 
The  contact  time  is  10.6  m  sec,  and  the  trace  is  taken  on 


Reproduced  from 
best  available  copy. 


%  TRANSMISSION 


The  linear  velocity  was  measured  using  an  oscilloscope 
(Tectionics  type  556)  to  monitor  the  pulsed  reference 
signal  from  the  microwave  generator.  The  scope  also  monitored 
the  NO  +  0  glow  AC  signal  at  a  fixed  port  when  NO  was  put 
in  the  system.  Since  the  production  of  oxygen  atoms  was 
pulsed,  and  for  these  observations,  the  duty  cycle  was 
reduced  to  10%,  the  NO  +  0  glow  signal  was  also  pulsed.  The 
delay  between  the  production  of  the  oxygen  atoms  and  the 
observation  of  the  signal,  which  could  be  viewed  as  the  delay 
between  the  reference  and  the  glow  signal,  was  observed  on 
the  scope.  This  delay  was  related  to  the  velocity  in  the 
reactor.  Measurements  at  two  different  ports  would  give 
the  delay  time  between  JUlses,  and  hence  the  velocity.  A 

typical  velocity  vs.  distance  profile  for  the  system  is  given 
in  Figure  4. 
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Ill  I  1 

10  20  30  40  50  60 

Distance  (cm) 

Figure  4.  Velocity  vs.  distance  downstream  from  the  point  of 
reactant  mixing.  The  total  pressure  is  0.96  Torr.  The  error 
in  the  measurement  is  greatest  between  the  first  two  windows. 
Thus,  the  velocity  appears  to  be  essentially  constant  down 
the  flow  tube,  which  implies  a  constant  density  system. 

The  velocity  measurement  made  is  the  centerline  velocity  from 
bolus  washout  theory.  The  average  velocity  is  half  of  this 
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value  . 


THE  MECHANISM 


Gene  ral 

The  reaction  of  oxygen  atoms  with  ethylene  has  been 
phenomenologically  found  to  produce  infrared  radiation,  at 
least  partly  in  the  2.7y  region.  The  major  effort  of  this 
work  is  devoted  to  the  quantitative  explanation  of  that  radi¬ 
ation  in  forms  of  the  detailed  molecular  mechanism.  The  gen¬ 
eral  ground  state  system  of  intermediates  resulting  from  the 
combination  of  0  and  C^4  weH  established  primarily  from 
the  discharge-flow-tube-time-of-flight  mass  spectrometric  study 
of  Niki  et  al  (1969)  and  the  direct  observation  of  the  products 
CH3  and  CHO  in  the  high  intensity  crossed  beam  experiment  of 
Kano f Sky  et  al  (1972).  Th*s«  two  studies  along  with  the  simi¬ 
lar  flow  tube  mass  spectrometric  work  of  Herron  and  Penzhorn 
(1969)  have  established  that  the  primary  elementary  step  (901) 
in  the  reaction  mechanism  is 

0  +  C2H.  +  CH3  +  CHO  (1) 

Other  minor  species  were  observed  by  Kanofsky,  et  al  (1972) , 
but  notably,  the  production  of  CH2  and  H2CO  were  not  observed 
to  occur  in  the  primary  step.  The  eventual  stable  products 
formed  by  secondary  reactions  include  H2CO,  CO,  H2,  02 ,  and  H. 
Consideration  of  the  temporal  behavior  of  these  products  as  well 


* 


Numbering  refers  to  Table  II. 
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as  isotopic  experiments  led  Niki,  et  al  (1969)  to  postulate 
that  the  following  reactions  occur: 


CHj  +  0  -*■  HCHO  +  H 

C2) 

CHO  +  0  -**  CO  +  OH 

(4) 

HCHO  +  0  -*•  CHO  +  OH 

(17) 

CHO  +  0  +  C02  +  H 

(5) 

OH  +  0  +  02  +  H 

(6) 

CHO  +  H  -*■  H2  +  CO 

(7) 

Water  was  not  observed  to  be  a  product  of  this  system  by  Niki, 
but  his  sy;. -  cm  employed  very  small  concentrations  of  C2H4  with 
0  in  excess,  so  that  the  OH  concentration  and  hence  toe  rate 
of  the  reaction: 

OH  +  OH  -►  h2o  +  0  rn) 

were  both  small.  In  a  more  general  case,  H^O  might  be  a  sig¬ 
nificant  product. 

Both  OH  and  H20  are  produced  by  exothermic  reactions  in 
the  above  sequence  under  conditions  where  each  is  the  new  bond 
formed  by  a  neutral  rearrangement  reaction.  Under  conditions 
such  as  this,  it  is  well  established  that  a  significant  frac¬ 
tion  of  the  exothermicity  of  the  reaction  can  appear  as  vib¬ 
rational  energy  in  the  new  bond  (Smith  (1972);  Pacey  and  Polanyi 
C 19 7 1 )  ;  McGrath  and  Norrish  (1957);  Parker  and  Pimentel  (1969); 
Jonathan,  et  al  (1971);  Stair  and  Kennealy  (1967);  Hushfar,  et 
al  (1971)  and  many  others). 
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TABLE  II 


Reaction 

Rate 

Constant , 
300  °K 

Reference 

AU  k  cal 
r  gmole 

1 

C,H4  +  0  -►  CH3  +  CHO 

8 . 42 " 13 

Davis,  et  al 
(1972) 

-42.45 

2 

CH3  +  0  HCHO  +  H 

3 . 0  " 11 

Niki,  et  al 
(1969) 

-66.54 

3 

HCHO  +  0  ->  CHO  +  OH^ 

Calculated 

4 

CHO  +  0  -*■  CO  +  OH 

l.O'12 

Niki,  est. 

-  73.09 

5 

CHO  +  0  -»•  C02  +  H 

3.0'13 

Niki,  est. 

-98.05 

6 

OH  +  0  -»■  02  +  H 

5. O'11 

Kaufman  (1964) 

-16. 79 

7 

CHO  +  H  -*■  CO  +  H2 

4. O’12 

Niki,  et  al 
(1969) 

-75.62 

8 

°H=f(v=i)  -  OH  +  hv 

3.3 

Potter,  et 
al  (1971) 

9 

CHO  +  02  -»■  CO  +  H02 

2x10 ’ 13 

-18.88 

10 

OH^  +  0  +  02  +  H 

5.0’11 

Kaufman  (1964) 

11 

OH  +  OH  -+  H20  +  0 

2 . 5  "  12 

Dixon -Lewis  et 
al  (1966) 

-16  95 

12 

OH  +  OH  +  H20^  +  0 

2 . 5 ' 1  3 

Estimate 

13 

14 

CH3  +  02  +  M  -  CH302  +  M 

H  +  02  +  M  -*■  H02  +  M 

8 . 0 "  32 

3. O'32 

Heicklen 

Kaufman 

(1964) 

-47.4 

15 

H20^  +  M  -*■  H20  +  M 

l.O'12 

Estimate 

16 

H20^  -►  H20  +  hv 

30 

Penner  (1959) 

17 

HCHO  +  0  -*■  CHO  +  OH 

1 . 6 " 1 3 

Niki ,  et  al 
(1969) 

-24.77 

18 

CHO  +  0  CO  +  OH^ 

Calculated 
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OH  <  HCHO  -*•  H20'  +  CHO 
OH  +  H2  •+  H20  +  H 
OH  +  H02  -►  H20  +  02 
H02  +  0  OH  +  02 
H  +  H02  OH*  +  OH 
OH  +  H02  -*■  H20^  +  02 
C2H4  +  0  -  CH3  +  CHO* 
CHO*  +  0  -v  CO  +  OH* 

H  +  HOr*-  OH  +  OH 
0  +  H02  •>  OH*  +  02 
CO  +  OH  +  C02  +  H 
OH^  +  02  -*•  02  Cb 1 E  )  OH 
OH*  +  H  ■>  h,  +  0 
02  (b  E)  +  02  -*■  0,  +  02 
02  (b1!!)  +  0  -  o2  +  0 

°2  (b1!)  +  H  -*■  02  +  H 
H  +  HCHO  H2  +  CHO 

CHO*  +  02  -*■  CO  +  H02 
0  ->  1/2  02 
H  -►  1/2  H2 
OH*  •>  OH 
h2°^  H20 

CHO*  -*  CHO 

°2  (blV  °2  (*3Z) 

Dummy 


1.4'12 

Estimate 

6.6'15 

Greiner  (1969) 

i.o-11 

Kaufman  (1964) 

1 .o"11 

Lloyd  (1970) 

1 .0 " 12 

Crutzen  (1971) 

l.O'12 

Estimate 

5.2xl0'14 

Estimate 

l.O'11 

Estimate 

• 

o 

1 

h-» 

h-» 

Estimate 

l.O’12 

Estimate 

1.5-13 

Wilson  (1972) 

3x10 " 14 

Potter,  et  al 

lxl. O’12 

Estimate 

3xl0"16 

Estimate 

2xl0"13 

Estimate 

2x10  13 

Estimate 

1x10 " 13 

Westenbure  and 
Dehaas  (1972) 

1 xlO  11 

Estimate 

137 

Estimate,  wall 
deacti va  tion 

137 

Estimate 

137 

Estimate 

137 

Estimate 

137 

Estimate 

137 

Estimate 

-119.234 

-  71.87 

-  54.2 
-38.4 


OH*  +  M  -*•  OH  +  M 

OH*  +  H2  ->  H20  +  H 

1 . 5x10  14 

Estimate 

lx] 0 " 1 4 

Estimate 

C2H4  +  OH  -*>  CH3  +  HCHO 

OH*  +  C2H4  + 

1 . 8x10' 12 

Morris  6  Niki  (1971) 

8- 


Thus,  in  the  general  case,  both  the  Oil  fundamental  and  the 
H20  v3  band  could  be  expected  to  contribute  to  emission  near 
2 . 7y .  In  order  to  simplify  the  .nitial  kinetic  system,  it 
was  decided  to  work  at  small  C2Hj  c°ncentration  in  order  to 
minimize  the  H20^  production.  In  addition,  the  02  concentra¬ 
tion  was  minimized  in  order  to  avoid  the  reactions 


H  +  02  ♦  M  *►  H02  +  M 

0  +  H02  -*•  OH^  +  02 

H  +  H02  -*•  OH^  +  OH 

OH  +  H02  -  H20+  +  02 
■ 

The  minimization  of  the  C-H. 

2  4 


(14) 

(29) 

(24) 

(25) 

also  prevented  significant  con- 


liiuuliuns  fxum 


(47) 


OH  +  C2H4  -*■  products 
OH  +  H2  -►  H20^  +  H, 

as  well  as  minimizing  the  contribution  of 
OH^  +  C2H4  -*■  products 

0HT  +  HC  -  fragment  -*■  products 

A  , 

to  the  OH  loss.  Under  such  circumstances,  the  major  OH*  loss 

process  will  be  via 


(48) 


OH 


^  +  0  -*•  07  +  H, 


(10) 


which  is  extremely  fast,  even  for  the  ground  state  OH  molecule. 
In  fact,  reaction  w'ith  0  as  a  homogeneous  loss  process  i  so 
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rapid  that  0HT  deactivation  on  the  wall,  which  would  be  diffusion 
controlled  even  if  every  wall  collision  were  effective  in  de¬ 
activation,  can  be  neglected.  However,  it  was  found  that,  in 
the  general  case  small  additional  loss  mechanisms  provided  by 


OH^  +  02  -*•  02(1Eg)  +  OH  (31) 

and 

OH^  +  M  -*■  OH  +  M 

were  needed  to  fit  the  observed  data  at  long  contact  times. 
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The  problem  of  analyzing  the  relative  role  of  each 
reaction  in  the  complete  mechanism  of  Table  II  was  approached 
both  numerically  and  analytically.  In  the  tormer  method,  all 
the  reactions  were  included  and  the  resulting  differential 
equations  were  integrated  using  the  Keneshea  (1967)  code  which 
incorporates  the  Runge  Kutta  Merson  integration  procedure  with 
a  steady  state  linearization-iteration  algorithm  which  is 
implemented  when  the  equations  become  stiff.  This  is  a  standard 
procedure  previously  employed  for  such  problems  (Kummler  and 
Bortner,  196«;  Gutowski ,  1971;  Roginsky,  1972).  The  numeri¬ 
cal  integration  was  employed  for  two  purposes.  First,  Hypotheses 
regarding  the  relative  importance  of  the  various  reactions  and 
species  could  be  readily  tested  under  prototype  conditions  and 
the  experimental  conditions  could  be  adjusted  to  optimize  the 
desired  measurement.  Second,  consistency  between  the  final 
predicted  and  calculated  intensities  could  be  obtained, 

It  was  found  that  H20^  emission  strongly  increased 
when  the  C2H ^  concentration  in  the  system  exceeded  a  few  microns, 
and  became  insignificant  at  concentrations  of  C2H^  less  than  ly 
as  shown  in  Figures  5  and  6  from  the  numerical  model.  Using 
the  capillary  tubing  provided  the  needed  restriction  to  the 
C2H^  flow  into  the  reactor  and  allowed  the  introduction  of 
C2H4  partial  pressures  as  low  as  O.ly  at  a  1  tori  total  pressure 
in  the  flow  system. 
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figure  5.  The  log  of  the  concentrations  of  the  species  of 
major  interest  at  high  C2H4  concentration  as  a  function  of 
reaction  time. 


TIME, SEC  (X101  ) 


Figure  6.  The  major  species  concentrations  at  low  C2H 
concentration  as  a  function  of  reaction  time. 


Assuming  experimental  conditions  preventing  the  forma¬ 
tion  of  H20*  and  favoring  OH*  as  in  Figure  6,  we  can  examine 
the  reaction  mechanism  as  shown  in  Figure  7.  Oxygen  atoms  plus 
ethylene  yield  the  methyl  and  formyl  radicals,  CH-  and  CHO. 

4  3 

Then,  0HT  production  occurs  via  two  routes.  The  first  is  the 
immediate  reaction  of  the  formyl  radical  produced  in  the  pri¬ 


mary  step: 


0  +  CHO  CO  +  OH^ , 


(18) 


which  is  extremely  rapid  (Browne,  et  al,  1971).  CHO  is  also 
lost  by  the  rapid  process: 


CHO  +  0-,  products 


(19) 


which  will  dominate  at  high  02  concentrations.  With  the  steady 
state  postulate  on  CHO,  for  purposes  of  this  data  analysis,  the 
rapidity  of  this  reaction  uermits  the  reaction  to  be  conceptu¬ 
alized  in  the  following  way.  The  reaction  mechanism  is  equiva¬ 
lent  to  the  introduction  of  CHO  rather  than  C^  since  the 
steady  state  immediately  converts  the  CHO  to  OH*  with  no  time 

lag.  The  rapid  steady  state  also  occurs  for  the  methyl  radical 
which  is  lost  via 


0  +  CH3  +  H2CO  +  H 


The  result  of  the  oxygen  atom  plus  methyl  radical  does  not, 
however,  result  in  direct  production  of  OH*.  As  shown 
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in  Figure  7,  the  OH^  results  from  the  relatively  slow 
subsequent  process 

0  +  H2C0  +  CHO  +  OH^,  (  3) 

and  the  CHO  also  reacts  with  0  as  above  to  produce  an  addi¬ 
tional  OH^.  Because  reaction  3  is  slow,  all  OH^  produced  via 
the  CHj  primary  product  will  exhibit  a  time  lag,  whereas  the 
OH^  produced  from  the  formyl  radical  which  is  a  primary  pro¬ 
duct  will  exhibit  no  apparent  time  lag-  This  suggests  a  method  of 
obtaining  both  quantum  yields  from  analysis  of  the  radiation 
as  a  function  of  flow  tube  position.  While  the  numerical  solu¬ 
tion  illustrated  in  either  Figure  5  or  6  indicates  a  rapid  rise 
of  0HT  and  CHO  from  zero  at  zero  time,  extrapolation  of  the  Oil' 
curve  indicated  by  the  dotted  line  in  Figure  6  will  yield  an 
intercept  which  represents  the  01 1^  production  due  to  the  hypo¬ 
thetical  introduction  of  CHO  at  zero  time.  This  intercept,  as 
shown  below,  permits  the  obtention  of  the  rate  constant,  k^g, 
for  the  reaction 

l 

0  +  CHO  -*■  CO  +  0HT.  (18) 

Analysis  of  the  data  for  long  contact  times  permits  estimation 
of  the  rate  reaction  3, 

0  +  CH20  -V  OH^  +  Clio.  (  3) 

This  division  of  time  zones  is  supported  by  the  work  of  Kriegcr, 
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0XV6FN  ATOM  ATTACK  OM  £THYl£rt£ 


Figure  7.  Schematic  of  the  Reaction  Mechanism  Proposed 


Malki  and  Kummler  (1972) ,  who  found  that  the  OH  Meinel  bands 

I 

(v  ^  7  to  9)  were  produced  early  in  the  flow  system  as  shown 
in  Figure  8. 
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short  contact  time  O  1-3  m  sec)  taken  from  Krieger,  Malki  and 
Kummler  (1972).  The  total  pressure  is  1.28  Torr  and  the 
hydrocarbon  partial  pressure  is  20  millitorr. 
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First,  the  steady  state  assumption  will  be  invoked  for 

A  A 

OHT  because  the  time  constant  for  the  loss  mechanism  for  0HT, 

OH^  +  0  k^°  02  +  H,  CIO) 

is  small  compared  to  the  reactor  contact  time,  k^g  may  be 
estimated  in  two  ways.  First,  we  will  assume  that  the  lower 
bound  is  the  rate  constant  for 

OH  +  0  ->•  CL  +  H  (6) 

as  determined  by  Kaufman  (1964)  and  by  Clyne  (1963) :  k^g  = 

5  x  10  ^  cc/sec.  Thus,  at  a  typical  oxygen  atom  number  density 
of  ?  X  1014  1/rC;  thf»  pffprl-i  VP  fi.rqt  nrrlpr  r?rp  rnn^^t  for 

4  4  4 

the  loss  of  OH r  is  10  ,  or  the  time  constant  for  OHT  loss  is 

0.1  msec,  compared  to  the  contact  time  at  the  first  window  of 
2,42  msec.  Other  potential  loss  mechanisms,  such  as 

0H+  +  OH  -*•  H20  +  0 

and 

OH^  +  H2  *»■  H20  +  H 

OfP  +  OH  +  hv  (8) 


are  at  least  10  times  slower  under  most  experimental  conditions 
employed  in  the  present  study. 


The  upper  limit  for  k^g  would  be  the  gas  kinetic  rate 


constant, 
k  - 


where  for  this  system 


°OH  "  3A’  M0H  =  17 

aQ  ■  1.3A,  *  72,  R  «  8,3  x  10 8 ,  T  =  300°K 

*  *  *  a0H-0  “  2*2A»  M  =  . 2 36 " 1 
*10  -  1.7  x  10'10  cc/sec 

The  lower  limit  is  taken  to  be  the  value  of  k^g. 

Thus,  from  Table  II,  in  the  steady  state, 

i  VH2C0]  [°]+ki8£CH°)  [0]+k24  IH]  [H02]+k29[0]  [H02] 

OH  -  - 1c10  10  J -  (1) 

Reaction  24  is  found  not  to  be  important  in  the 

mechanism;  An  experiment  that  prcved  the  ineffectiveness 
of  reaction  29  was  carried  out.  Molecular  ground  state 

oxygen  that  would  enhance  the  production  of  HO-  via  reac- 
tions  9  and  14  did  not  increase  the  0HT  (2.7p)  emission  as  reac¬ 
tion  29  suggested.  One  could  therefore  conclude  the  quantum 
yield  of  reaction  23  was  not  appreciably  important.  Thus,  equa¬ 
tion  (1)  reduces  to: 

i  k,  •  CH-0  +  k i  -  *  CilO 

OH*  «  - - — rr - ^ -  (la) 

K10 
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Under  these  assumptions  and  the  postulate:  (kg)  .  0 2  >  (k^  +  k^) 
which  is  the  only  condition  consistent  with  the  data,  we  can 
compute  the  OH^  concentration. 

Neither  k4 ,  kg,  nor  kg  are  well  known.  The  above 

reduces  the  steady  state  equation  for  CHO  to 


k,  C,H.  4  k,  7  Crl.,0 
rH0  .  (-1  2  ^  •  0 


(8) 


The  molecular  oxygen  concentration  is  constant  down  the 
tube  as  described  by  the  computer  model.  This  allows  us  to 
write,  upon  substituting  the  time  dependency  in  equation  (8): 


CHO 


,  .  u  ,  r  -kj/Odt  -k, 7/0dt 

=  klC2H4lp{e _ «(-2k17+k1)4:c},/e  }  •  0  (9) 

k  1  kg  Q 


where 


k  1  =  kl  "  k17 


Grided  by  the  numerical  integration  that  describes  this 

approach  the  formaldehyde  concentration  is  considered  to  be  con¬ 
stant  down  the  tube.  It  is  reasonable  to  predict  a  small, 

slowly  increasing,  contribution  to  the  emission  from  this  species 
according  to 


I 


due  to 

ch2o 


kl C2H4  I  0  ,  -•'rk17,0-dt 

i  -  ( e 


- / k. '0' dt 

e  ) 


A  k- 
m  3 


(10) 


Ibc  exponential  function  of  the  second  term  is  larger  than  that  of 

of  the  first.  This  should  give  a  small  upward  trend  in  the 
intensity . 
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'  This  small  rise  is  not  significant  when  compared  to  the 
contribution  from  CI10 


K . 'k..  Q *k- *C-H  .  [ n  -k. /Odt 

*  I  due  to  1  ■  .  o  <V»17> 

CHO  k10’k  lk9°2 


-k17/Odt  -k  t 

e  17  •  k1?3  •  0|oe  * 


where  k^  =  wall  recombination  (of  oxygen  atoms')  rate  constant. 


Both  the  increase  of  /Odt  and  the  decrease  of  the  oxygen  atom 
concentration  due  to  wall  recombination  reduces  the  intensity 
as  a  function  of  time. 

This  is  best  represented  by  the  expression 


iogjc  I  =  log10  A  -  kwt 

where 


(11a) 


The  time  dependency  of  A  i$  small  ard  will  give  an  approximate 
log  linear  equation  over  the  range  of  interest.  This  explains 
the  semi  log  straight  line  relationship  which  appears  to  fit 
experimental  data. 
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The  approximation  of  constant  / Odt  does  not  effect  the 
accuracy  of  the  calculated  k.g.  It  is  illustrated  earlier 
that  an  intercept  when  t  =  0  is  a  valid  concept  and  is  equiv¬ 
alent  to  a  hypothetical  introduction  of  CHO  at  time  zero. 
Equation  12  reduces  to 


log1Q  I  =  log 


10 


Ank18klC2H4 I 
”k10*k9°2 


0 


and 


A  k18  _  1  k10  °2 

n  •  E7  -  O'k,  up;;!-  (13) 

At  the  end  of  the  tube  when  CHO  production  decreased  to 
a  great  extent  due  to  the  wall  recombination  of  the  oxygen 
atoms,  the  only  significant  radiation  would  be  due  to  formal¬ 
dehyde.  This  argument  could  be  used  to  place  an  upper  limit 
on  Kj  The  intensity  at  the  last  window  is  taken  to  be 

the  limit  corresponding  to  the  highest  possible  k^  (from  equa¬ 
tion  11)  : 


k 


3 


I 


kl C2FI4  I  0 


limit  k  1 
- /k.  70Jt 
(e  17 


k 


e 


10 

-  7k"0dt 


where 


(14) 


1 1  limit  =  the  maximum  intensity  due  to  CIIo0 

A  plot  of  intensity  vs.  C7H4  at  constant  contact  time  is 
obviously  expected  to  give  a  straight  line  relationship.  The 
curvature  in  the  experimental  data  presented  in  the  results 
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section  is  attributed  to  the  small  dependence  of  the  oxygen 
atoms  on  ethylene.  This  dependence  becomes  least  important 
at  low  ethylene  concentration.  Thus  an  initial  slope  of  the 
data  is  a  good  approximation  of  the  slope  described  by  equa¬ 
tion  c  11)  : 


Slope 


An*kia*ki  -ki/Odt 

18  1  ■  OL  (e  1  (k,  -  2k17)  ♦ 


l'k9°2 


0 


-k1?/Odt  -k  t 

e  •  k1?j  •  e 


(15) 


The  term  / Odt  is  calculated  numerically,  k  is  the  wall  recom- 
.  w 

bination  rate  constant  of  oxygen  atoms.  kw  is  calculated  by 
observing  the  NO  +  0  glow  signal  decay  down  the  tube. 
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Results  and  Discussiun 


It  nas  been  proposed  that  the  reaction 

CHO  +  0  -*■  CO  +  OH^  (2 . 7p)  (-18) 

is  primarily  responsible  for  the  emission  in  the  flow  tube  sys¬ 
tem.  It  has  also  been  postulated  that  the  reaction 

CH20  +  0  -*■  CHO  +  OH^  (7.7)  (  3) 

would  predominate  at  longer  contact  times  due  to  the  decrease 

of  the  oxygen  atom  concentration  (see  equation  10  and  11  in  the 
previous  section). 

Under  this  mechanism,  in  which  CHO  is  primari1''  ^s1,  !y 
reaction  with  02 ,  it  is  predicted,  as  shown  in  Figures  9  and  10 
that  the  observed  intensity  should  be  nearly  proportional  to  the 
C2H4  concentration,  until  C.,!^  reacts  significantly  with  and 
therefore  depletes  the  0  atom  concentration. 

In  Figures  9  and  10,  the  solid  line  represents  the  theoreti¬ 
cal  prediction  of  the  computer  model  and  the  points  represent  the 
experimentally  obtained  OH*  concentration  using  the  actinometric 
standard  of  the  NO  +  0  glow  as  determined  by  Fontijn,  et  al.  The 
rate  constants  for  the  computer  fit  are  obtained  using  the  analytic 
solution  for  Postulate  3,  based  on  the  data  of  June  15  (Figure  9) 
as  previously  described.  The  computer  code  is  then  used  to  pre¬ 
dict  the  data  of  June  12  which  has  a  higher  oxygen  atom  concentra¬ 
tion  as  well  as  a  higher  02  concentration  (Figure  10). 
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3.91  x  10AH/cc 


Contact  Time 
2.42  m  sec 


3.58  m  sec 


6.00  m  sec 


8.28  m  sec 
10.54  m  sec 


16.6  m  sec 


The  experimental  data  shows  some  deviation  at  higher 
ethylene  concentrations  and  short  times  where  significant  devia¬ 
tions  from  the  proposed  mechanism  are  expected  to  appear  because 
of  H20*  emission,  but  the  data  generally  verify  the  expected 
proportionality  between  observed  intensity  and  initial  C2»4  con¬ 
centration  over  an  order  of  magnitude  in  concentration  and  time. 

A  plot  of  (intensity)  vs.  time  should  enable  two 

pieces  of  information  to  be  obtained:  a  value  for  and  an 
upper  limit  on  k,.  Figure  11  represents  the  typical  plot  of  the  inten- 
sity  of  each  point  divided  by  the  corresponding  C2H4  values  as  a 
function  of  time,  since  the  equations  predit  essentially  a  first  order 
dependence  on  C^.  The  slopes  in  the  plots  primarily  describe 
the  wall  recombination  rate  constant  and  are  weakly  dependent  on 
C2H4  and  /0dt  ( see  equation  11  of  the  previous  section). 

The  time  interpretation  of  the  intercepts  of  Figures  10-11. 
has  justified  plotting  equation  (11a)  in  the  described  fashion. 


Thus 


where 


r't=( 


AJt  =  0 

W7 


^t-0  = 


1 1  =  0 


An‘k18'kl 
^10  *k9  *°2 


IU-0  1S  the  ratio  of  intensity  to  concentrati 


on 
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A|j._q  is  defined  in  equation  12,  and  is  not  timp  dependent 


k18  a  1  1 t=0  *kr/°2 

k  9  0 • *An 

where  A  is  the  Einstein  coefficient  associated  with  reaction  18. 
n  k 

o 

The  obtained  value  for  r —  =1.7  supports  the  earlier  postulate 

9  k18 

that  kg  C>2  >>  kjg*0,  or  0^  >>  ^ —  •  0  where  the  atomic  oxygen 
concentration  is  8p  in  80y  of  0^  . 

From  the  maximum  Cl^O  intens ity  observed  k^  was  found  to 
be  3  x.  10  .  The  intensity  at  the  6th  window  was  taken  to  be 

the  highest  that  reaction  (3)  could  produce^ 
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a^d  Am  is  the  associated  Einstein  coefficient  with  reaction  3. 

The  step  calculations  done  for  the  above  results  are 
represented  in  Table  IV. 

The  values  for  A  and  used  in  this  work  were  assumed 

to  be  10  1/sec  and  3.3  1/sec  for  reactions  18  and  3,  respec¬ 
tively,  guided  by  the  exothermicitv  of  the  respective  reactions. 
It  is  possible  that  the  OH ^  emission  is  coming  from  a  transition 
higher  then  the  (1-0)  transition,  via  reaction  18,  while  the 
smaller  exothermicity  of  reaction  3  would  justify  using  3.3  1/sec 
(Potter,  et  al ,  1970)  . 
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Figure  11.  Log10  Intensity  as  a  function  of  time.  Bach 
point  was  divided  by  its  corresponding  ethylene  concentration t 
upon  doing  this  almost  all  points  corresponding  to  the  same  time 
overlapped.  The  slope  of  this  curve  mainly  described  the  deple* 
tion  of  oxygen  atoms. 
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Thus,  in  summation,  we  have  observed  2.7u  radiation 
emunating  from  the  0  +  reaction  system  which  we  attribute 

primarily  to  the  reaction 

DIO  +  0  OH^  +  CO  fl8) 

OH^  is  lost  principally  by  reaction  with  oxygen  atoms,  through 
OH^  +  0  +  02  +  H;  (10) 

CHO  must  be  lost  principally  via 


CHO  +  02  products 


(  9) 


None  of  the  rate  constants  needed  here  can  be  obtained  uniquely 

from  this  study.  However,  the  critical  parameter  in  the  studv: 

)  »  can  be  determined  by  matching  the  experimental  data 
10  9  ]< 

extrapolated  to  zero  contact  time  with  the  value  .  .  Assuming 

1  k10  k9 

that  OH  reacts  at  the  same  rate  as  OH  in  the  essentially  gas 

kinetic  reaction  (10)  producing  02 ,  we  find  that  klg/k9  =  1.7  f  0.2 

Using  this  ratio  and  the  error  limits,  we  can  reproduce  all  the 

data  theoretically  to  within  the  error  limits. 
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Part  II:  THEORETICAL  DETERMINATION  OF  VIBRATIONAL  EXCITATION 

OF  MOLECULES 


R.  Marriott 


aEORETICAL__  DETERMINATION  OF  VTPPftTTONAX,  EXCITATION  HP  ang 
The  theoretical  program  at  RIES  has  been  primarily  con¬ 
cerned  with  the  calculation  of  cross  sections  for  the  vibra¬ 
tional  excitation  of  molecules  by  heavy  body  collisions. 

This  is  carried  out  using  a  numerical  close  coupling  solution 
of  the  wave  equation  for  the  colliding  system  based  on  the 
breathing  sphere  model  of  Schwartz,  Slawsky,  and  Herzfeld  (1954) 
to  represent  the  target  molecular  vibrator  and  an  empirical 
spherically  symmetric  scattering  potential,  but  neglecting 
any  structure  in  the  impacting  particle. 

For  many  gas  mixtures  of  interest  direct  vibration  - 
vibration  (W)  energy  transfer  processes  are  also  important. 
Nevertheless  no  calculations  of  such  processes  to  the 
precision  afforded  by  the  close  coupling  approximation  have 
been  published  up  to  the  present  time. 

The  RIES  molecular  collision  code  has  now  been  extended 
to  include  vibrational  states  in  the  impacting  particle. 

This  in  turn  permits  both  the  calculation  of  significant  W 
cross  sections  simultaneously  with  the  VT  cross  sections 
within  the  same  breathing  sphere  approximation,  and  the 
estimation  of  the  effect  produced  on  target  particle  VT 
cross  sections  by  coupling  in  the  impact  particle  states. 

Calculations  on  the  CO  -  CO  system  are  currently 
underway  to  study  the  magnitude  of  VV  cross  sections  for 
states  in  exact  energy  resonance  and  to  determine  the 
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effect  of  the  impact  particle  states  on  the  calculated  VT 
cross  section  for  the  excitation  of  the  first  vibrational 
level  of  CO  from  the  ground  state.  This  cross  section 
dominates  the  vibrational  relaxation  of  CO ,  and  since  the 
previous  VT  results  were  in  good  agreement  with  experiment, 
the  relevant  cross  section  would  not  be  expected  to  be 
particularly  sensitive  to  the  addition  of  impact  particle 
states . 

Preliminary  results  for  collision  systems  in  which  coupling 
of  both  2  and  3  states  in  the  impacting  and  target  particles 
has  been  retained  are  shown  in  Table  I. 

The  cross  section  for  the  excitation  of  the  first  vibra¬ 
tional  level  is  reduced  by  a  factor  of  2.5  from  the  correspon¬ 
ding  3  state  calculation  which  ignores  the  structure  in  one 
purticle .  The  VT  cross  section  is  reduced  by  a  factor  of  4  from 
the  original  4  state  results  which  give  agreement  with  vibra¬ 
tional  relaxation  data.  If  a  similar  reduction  occurs  at  all 
energies  it  would  significantly  affect  the  comparison  of  theory 
with  experiment.  However  it  can  be  seen  that  there  are  now 
additional  competing  process  which  must  be  taken  into  account. 

The  energy  exchange  cross  sections  are  very  large.  Clearly 
processes  of  the  (V-V)  type  C0(0)  +  C0(2)  C0(1)  +  C0(1) 

will  be  orders  of  magnitude  more  effective  in  depopulating 
the  higher  vibrational  states  of  CO  than  (V-T)  processes  such  as 
CO(O)  +  CO (2)  +  CO ( 0)  +  CO (1)  +  . 266eV. 

The  magnitude  of  the  V-V  cross  sections  in  comparision 
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with  V-T  collisions  is  found  to  be  caused  by  not  only  an 
increase  in  the  individual  partial  cross  sections,  but  also 
by  the  fact  that  of  the  order  of  twice  as  many  (300:150) 
partial  waves  contribute  significantly  to  the  (V-V)  total. 
This  re-emphasizes  the  long  range  nature  of  the  V-V  inter¬ 
action  in  a  very  concrete  way. 

While  the  results  obtained  so  far  seem  reasonable, 
they  must  be  regarded  as  tentative  since  it  is  evident  that 
a  sufficient  number  of  molecular  states  to  ensure  convergence 
of  the  cross  sections  lave  not  yet  been  coupled  together. 

It  was  shown  in  the  original  CO-CO  (V-T)  calculation 
that  4  or  5  states  had  to  be  included  for  convergence. 
Unfortunately,  as  presently  constructed,  the  molecular 
code  will  not  handle  the  case  of  virtual  states;  all  coupled 
states  must  be  energetically  accessible. 

This  means  that  for  any  specific  selection  of  excited 
states  th^re  is  a  certain  collision  energy  threshold  below 
which  the  code  cannot  be  run.  This  was  not  an  important 
restriction  for  the  V-T  calculation  where  only  one  excited 
particle  was  considered.  For  the  new  V-V  code  with  both 
particles  excited  however  it  is  more  serious.  Indeed 
including  5  states  in  each  CO  molecule  would  prevent  the 
code  being  run  at  a  collision  energy  less  than  2.7  eV, 
outside  the  region  of  interest  for  many  purposes. 

It  is  necessary  therefore  to  further  extend  the  code 


Table  I.  inelastic  cross  sections  (in  units  of  ira z)  for  CO-CO  collisions  at  a  total  system 
energy  of  1.5  eV.  ° 
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to  accept  the  case  of  virtual  states.  Although  this 
poses  no  analytic  problem,  indeed  close  coupling  electron 
scattering  codes  are  so  construe  Led.  it  will  require  a 
fairly  extensive  re-ordering  of  the  program  subroutines  and 
particularly  the  treatment  of  the  asymptotic  solutions. 

This  program  extension  is  now  underway  and  when  operational 
the  CO-CO  test  calculation  will  be  completed. 


A  molecular  beam  apparatus  donated  by  NASA  is  being 
used  to  measure  the  intermolecular  potentials  required  as 
input  for  close  coupled  calculations  of  excitation  cross 
sections.  A  number  of  improvements  to  the  apparatus  have 
been  necessary;  these  include  the  substitution  of  a  low 
energy  electron  bombardment  ion  source  and  improved  lens 
system  which  provides  a  stable  beam  with  very  low  energy 
spread.  In  addition  the  scattering  gas  handling  system  has 
been  reconstructed  to  eliminate  the  effects  of  surface 
outgassing  and  pressures  are  measured  with  a  capacitance 
manometer  calibrated  against  a  McLeod  gauge. 

Preliminary  experiments  with  the  Ar-Ar  system  yielded 
potentials  which  were  not  fully  in  agreement  either  with 
theory  or  with  previous  experimental  determinations.  Thus, 
in  order  to  check  that  the  apparatus  and  the  analysis 
programs  are  functioning  correctly,  we  have  measured  the 
potential  for  the  He-He  interaction,  for  which  highly  reliable 

I 

calculations  have  been  made  {1,2}.  For  this  system,  incomplete 
total  cross  sections  were  measured  as  a  function  of  relative 
energy  with  various  scattering  geometries. 

We  have  restricted  our  attention  to  the  case  where  the 
diameter  of  the  detector  p  is  larger  than  that  of  the  beam 
d  since  this  provides  more  signal  and  also  simplifies  the 
analysis.  Jordan  and  Amdur  have  shown  {3}  that  the 
measured  cross  section  is  a  maximum  when  d=p  and  decreases 
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asymptotically  as  d/p  tends  to  zero  or  infinity;  a  comparison 
of  curves  I  and  III  in  Figure  1  confirms  their  finding. 

These  were  measured  with  a  thermopile  ,  the  closed  symbols 
denoting  its  conventional  application  in  which  the  beam  transfers 
its  kinetic  energy  to  the  detector  and  produces  a  potential 
difference  across  the  junction  due  to  the  Peltier  effect. 

Such  measurements  require  proper  matching  of  the  low  source 
impedance  with  the  input  of  the  amplifier;  if  the  signal  is 
placed  across  a  high  impedance  input  then  the  secondary  electron 
emission  is  measured  (open  symbols  in  Fig.  1)  with  a  gain 
increased  by  a  factor  of  3000.  To  further  entend  the  range  of 
measurement  at  low  energies,  a  channeltron  electron  multiplier 
has  been  used  and  the  results  are  shown  as  curve  II,  Figure  1 
for  an  intermediate  value  of  d/p.  Inversion  of  these  curves 
pi o duces  potentials  which  are  virtually  indistinguishable 
from  the  best  theoretical  result  of  Yilbert  and  Wahl  {1},  and 
to  avoid  confusion  only  one  is  plotted  in  Figure  2.  An 
earlier  theoretical  result  of  Phillipson  {2}  is  similarly  close 

whereas  the  best  previous  experimental  results  of  Jordan  and 
Andrew  (3),  and  Kammer  and  Leonas  {4}  also  shown  in  Figure 
2  are  distinctly  lower. 


We  therefore  feel  confident  to  proceed  with  systems  for 
which  there  is  a  theoretical  check  and  are  currently  investi¬ 
gating  the  potential  for  the  Ar-Co  interaction.  It  is 
hoped  that  this  ,  together  with  the  potential  for  0  atoms  with 
exhaust  plume  species  will  soon  be  available. 
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